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Abstract — Resource allocation in orthogonal frequency division multi- 
ple access (OFDMA) networks plays an imperative role to guarantee the 
system performance. However, most of the known resource allocation 
schemes are focused on maximizing the local throughput of each cell 
individually, while ignoring the significant effect of inter-cell interference. 
Thus, the need of considering interference which in itself depends on the 
allocation of resources (i.e., subcarriers and powers) is evident. In this 
paper, the joint resource allocation problem via sum rate maximization in 
multi-cell uplink OFDMA is considered. The problem has a non-convex 
structure and is known to be NP hard even for single cell scenarios. 
Initially, upper and lower bounds are investigated. A centralized sub- 
optimal resource allocation scheme is then developed due to the inherent 
complexity of implementing the optimal solution. Furthermore, less 
complex centralized and semi-distributed schemes are proposed. The 
computational complexity of all schemes has been analyzed and the 
performance is compared through numerical simulations. Simulation 
results demonstrate that the semi-distributed scheme achieves comparable 
performance to the two centralized resource allocation schemes in various 
scenarios. 

I. Introduction 

In this paper, the problem of scheduling (i.e., determining the active 
users for a given time interval) and power allocation (i.e., assigning 
powers to these active users) via sum rate maximization in a multi- 
cell uplink OFDMA network is considered. Solving this problem 
in an optimal way is cumbersome due to the complexity issues 
involved in the joint optimization of a large quantity of continuous 
and discrete variables, non-convexity of the problem even for the 
continuous power variables, discrete nature of subcarrier allocations 
and the interdependency of power and subcarrier allocations. 

The sum rate maximization problem turned out to be more chal- 
lenging in the uplink scenario due to the distributed power constraints 
at each user. Simply allocating subcarrier to the users with best 
channel quality, the sum-rate may diminish as some active users may 
have considerably high channel gains but low transmission powers 
on a specific subcarrier than other competitors. Thus, for uplink, 
a subcarrier shall be assigned to the users who possess the largest 
transmission power-channel gain product in order to accomplish near 
optimal throughput in single cell scenario HI— 151. 

Recently, uplink OFDMA resource allocation via sum rate max- 
imization in single cell scenarios was investigated fll -t5l. In |T), 
a near optimal greedy scheme is proposed which is based on 
maximizing the marginal rate with iterative water-filling on the 
subcarriers. The authors in [2] generalized the strategy presented in 
171 by considering utility maximization and developed a polynomial 
time algorithm to compute an upper bound of the optimal solution. In 
(3j, sum rate maximization is studied with and without proportional 
fairness and a transmission power-SNR product (PSP) based ranking 
scheme is developed; however, the scheme is based on the statistical 
knowledge of the number of subcarriers assigned to each user. In (4), 
151 . the authors developed a sub-gradient based scheduling framework 
to compute the optimal solution of the relaxed problem. 

Some interesting papers in multi-cell uplink OFDMA networks 



are 161- 11101 . In J6), J7J, centralized and distributed schemes are 
studied which aims at minimizing the overall transmitted power 
subject to a rate constraint for each user. Nevertheless, most of 
the work is based on achieving a certain required rate (i.e., fixed 
rate adaptation) at each user rather than maximizing the network's 
throughput irrespective of any rate constraint that can exploit the 
spatial and multiuser diversity implied by the multicell dimensions 
in the OFDMA network. In (8), intercell coordination is used to 
maximize the network throughput by removing the cells that do 
not offer enough throughput; however, OFDMA networks need the 
use of multiple carriers to exploit the frequency diversity gain. In 
addition, some recent papers proposed low complexity distributed 
game theoretic solutions with certain pricing due to the apparent 
feasibility of implemeting them practically in the uplink OFDMA 
scenario |9), 1101 . Nontheless, all these approaches are suboptimal 
and no criterion is mentioned to caliberate their performance gap 
with respect to the optimal solution. 

Motivated by this fact, we investigate in this paper sum rate max- 
imization in multi-cell interference-limited uplink OFDMA network 
subject to an individual power constraint at each user. Initially, we 
investigate an upper (UB) and lower bound (LB) to the average 
network throughput. A centralized subcarrier allocation scheme is 
then developed assuming that a central controller exists possessing 
global knowledge of channel gains, interfering gains, and user powers 
in all the cells. Even though the centralized scheme has high 
computational and implementation cost, it can provide the guiding 
principles that help in designing low complexity centralized and semi- 
distributed schemes. In the semi-distributed scheme, we assume that 
the base stations (BSs) are allowed to cooperate in order to exchange 
interference information, subcarrier allocation decisions and jointly 
optimize power allocation in a distributed way. 

The rest of the paper is organized as follows: In Section II, the 
system model is defined and problem is formulated. In Section III, 
the bounds are derived and their complexity is analyzed. In Section 
IV and V, the centralized and semi-distributed resource allocation 
schemes are presented. Section VI demonstrates numerical results 
followed by concluding remarks in Section VII. 
Notation: Throughout the paper, we denote the sets of real and 
complex vectors of N elements by and C , respectively. 
Matrices are represented using boldface upper case letters while bold 
face lower case letters are used for vectors. AT(0, a 2 ) denotes a zero 
mean Guassian random variable with variance a 2 . 

II. System Model and Problem Formulation 

A network of L cells with a set of Ki users in each cell I is 
considered. Full reuse of the spectrum is assumed in all the cells, i.e., 
the reuse ratio is unity. Each BS is assumed to have N orthogonal 
subcarriers, and each subcarrier can be allocated to a single user per 
cell. The average throughput per cell is a function of both subcarrier 
allocation and power allocation variables. The sum rate maximization 



problem is formulated using the standard Shannon capacity formula 
G ni h t ,i = log 2 (l + 7n,fc,,i) where Cn,fc I> l,7n,fc I ,j represents the 
throughput and signal to interference plus noise ratio (SINR) of kj 
user at n th subcarrier in cell I, respectively: 

L K t N , h \ 

maximize a "- k ' lo S2 ( 1 + Vn 'Y. ) 0) 

" ' i=ik,=i«=i V ° + 1 "< l ' 

N 

Subject tO ^P„,fc,,i < Pfc,,max, Vfc( , VZ (2) 

71=1 

E^ a n ,ki,i — 1 ' Vn, VZ (3) 
k,=i 

a„, kl ,i £ [0,1], Vn,VZ,Vfc( (4) 

In (1), /„,; = Z)^=i,j/;Z)fc/=i a n,fc j ,iJ5n,fc j , 3 ffTi,fc J -,! which repre- 
sents the cumulative interference at n th subcarrier in cell I from 
users in all other cells, p n ,k t ,l denotes the power transmitted by fc* h 
user at the n th subcarrier in cell I, a n ,k t ,l represents the allocation 
of k\ h user at the n th subcarrier in cell I and h n ,k t ,l is the channel 
gain of kj h user at the n th subcarrier in cell I. Constraint (2) implies 
that the power spent by k\ h user on its allocated subcarriers cannot 
exceed the maximum power available denoted as max . For each 
cell, we pile the power allocation variables p n ,k l: l mto a vector 
Pn,i = [Pn,i,Pn,2, ■■■■,Pn,K l ] and then stack all the vectors into a 
power matrix P; of cell I where P; £ R JVxK ' . Constraint (3) restricts 
the allocation of a subcarrier to only one user. The channel gains 
hn,k t ,i and binary allocation variables a n ,ki,l are stacked up similarly 
in the matrices H; and A;, respectively, where A;,H; £ R JVxKi . 
Moreover, we define g n ,ki,lj as the interfering gain from the k\ 
user in cell I to cell j, Vj 7^ I at n th subcarrier. We pile these 
interfering gains into a vector g n ,y = [gr n ,i,« , 9n,2,ij 9n,K,,ij] 
and then stack all the vectors into a matrix Gy £ M. NxKl . In order 
to fully exploit the spatial and multiuser diversity gain offered by 
the multi-cell dimension, no quality of service (QoS) constraints are 
considered in (1). 

III. Bounds on the Network Throughput 

A. Lower Bound (LB) 

The LB for the optimum multi-cell throughput can be achieved 
when the inter-cell interference becomes maximum. Observing the 
dependency of inter-cell interference on the subcarrier allocation and 
power allocation variables, we assume each user in each cell is 
transmitting on each subcarrier with its maximum power. From (1), 
the average network throughput can be written as follows: 

1 L Kl N / h \ 

c(a,, p.) 4EEE I 1 + p ; fc 2 frf ) (5) 

The LB for the average network throughput taking the worst case 
inter-cell interference into account can be written as follows: 

c(a,,po > EEE^,' 1 * ( 1 + P ^h m, \ (6) 

!=U,=ln=l V / 

where ^ n>i = J2j=i j^i Z)fc.=i Pk j ,m.ax.g n ,k j ,l- The LB subcarrier 
allocations A; and power allocations P; in any cell / can be computed 
using Algorithm 1 where subcarriers are allocated to the users based 
on the Qi matrices where Qn,k t ,i = is a factor that 

takes into account the power-worst SINR product of each user at 



each subcarrier. These allocations are then utilized to determine the 
LB network throughput. Since, it has been shown that equal power 
allocation has insignificant performance loss compared to the optimal 
water-filling, power equalization among the allocated subcarriers is 
assumed |[T1. l3l. 

B. Upper Bound (UB) 

Establishing an upper bound on the uplink throughput in a multi- 
cell OFDMA network is significantly important in order to calibrate 
the performance of sub-optimal/near optimal resource allocation 
schemes. The UB can simply be computed by ignoring the effect 
of inter-cell interference in all the cells: 

c( Ai ,po <jEEE a »w% (i + p "' fc "ir fc "' ) (7) 

i = l fc;=l 71 = 1 ^ ' 

The UB subcarrier allocations A; and power allocations P; can be 
computed by substituting = 0 in Algorithm 1. These allocations 
can then be utilized to compute the average network throughput in 
(7). The near optimality of the average throughput is verified using 
the UB computed in [2] for the single cell scenario. However, the 
average network throughput as dictated by the UB allocations could 
be highly optimistic and may not be accurate for multi-cell systems. 
Consider a motivating example with two cells, two users and two 
subcarriers. Each user can transmit with a maximum power of 1W. 
Assume Hi= [1 0.9; 0.8 0.7] and H 2 = [1 0.9; 0.8 0.7]. Single 
cell allocation strategies that aim to maximize the local throughput 
of each cell suggest Ai,Pi and A2,P 2 =[1 0; 0 1]. Computing 
the UB using (7) results in 1.7655bps/Hz/cell. Now, assuming the 
knowledge of interfering link gains at each BS, i.e., Gi2= [0.9 0.2; 
0.2 0.9] and G 2 i= [0.7 0.1; 0.1 0.7]. Computing the throughput again 
while keeping the single cell allocations and taking into account 
the interfering gains leads to an average network throughput of 
1.1137bps/Hz. However, better allocations are possible if we consider 
Ai,Pi and A2,Pi = [0 1; 1 0] as per the criterion discussed in 
Section IV which enhances the resulting average network throughput 
to 1.5977bps/Hz. The throughput values calculated in the example 
may not be the maximum possible as the optimal powers are yet to 
be computed. This issue will be discussed in Section VI. 

Algorithm 1 Computing LB and UB Allocations in Cell / 

1) Input: [H;], [A;], [Pi], [Gy] where a„, fei ,; = 0, = 
[P fcl , m ax/iV] Vfc ; ,Vn 

2) For each fc* h user in cell I, power equalization is performed 
over all the subcarriers allocated previously to that user and 
the remaining unallocated subcarriers in the system. 

3) Using [Pi] from step 2, [H;] and [Gy], compute the matrix 
Q; for each cell /. 

4) Find the (n, fe;) pair that has the maximum value of Q n ,k lt i- 
Allocate subcarrier n to user fc;. 

5) Delete the n th subcarrier from the set of unallocated subcarri- 
ers. If there are still unallocated subcarriers in the system go 
to step 2, else terminate after distributing the maximum power 
at each user over all of its assigned subcarriers. 



Complexity Analysis: The (n,k{) pair at which the Q n ,ki,l In- 
comes maximum is allocated (Step 4), which has a complexity of 
a two dimensional search, i.e., O(KN). However, as soon as a 
subcarrier is assigned, each user updates its power as defined in 
Algorithm 1. This process iterates until all the subcarriers in all the 
cells are allocated and, thus, the time complexity of Algorithm 1 is 
0{KN 2 ). 



IV. Centralized Resource Allocation Schemes 

Assuming perfect knowledge of channel gains at a centralized 
controller, the optimal solution for (1) can be computed in the high 
SINR regime by an exhaustive search over all possible combinations 
of the allocations in each cell. For each possible allocation, optimum 
powers can be derived by transforming (1) into a geometric program 
(GP). Importantly, the power allocation problem is in itself a known 
non-convex problem for the general SINR regime. However, for the 
high SINR regime the problem turns into a convex GP problem. For a 
given set of allocation variables and considering a high SINR regime, 
the objective function in (1) can be rewritten as follows: 



maximize zl a ™M,^°& 

1 = 1 fc;=l n=l 



Pn,ki,lh n ,ki,l 
O 2 + IZ 



(8) 



Maximizing the SINRs is equivalent to minimizing the interference 
to signal ratio (ISNR), the objective can thus be rewritten as 



L K l N 

minimize ^2 z2 an M,^°& 

1 1 = 1 k; = l n=l 



3n,ki ,1 hn, kt ,1 



(9) 



Equivalently, (1) can be reformulated for high SINR regime as 
follows: 



Ki N 



V +In 



mimmize log 2 ]J II II I „ , , h 

1=1 ki=l ra=l v ' 1 ' ' 1 ' 

N 

subject to ^2p n ,k,,l < Pfc,,max, Vfc;,VZ 



(10) 



Note that the numerator in (10) is a posynomial and the denominator 
is a monomial, hence (10) is a GP problem in standard form that 
can be solved optimally through efficient interior point methods 
II II after performing the logarithmic transformation of variables 
1121 . However, even for small dimensions, it is not recommendable 
to compute the optimal solution due to the huge computational 
complexity of 0(K LN ) associated with an exhaustive search based 
subcarrier allocation phase. In addition, GP based power allocation 
method discussed above has two restrictions: high-SINR assumption 
and centralized time-consuming computations. These limitations, 
however, can be possibly reduced through the heuristics proposed 
in 1121 . Considering the high intricacy of implementing the optimal 
solution, a two-stage centralized scheme is developed. Moreover, a 
second centralized scheme with reduced complexity is discussed. 

A. Centralized Scheme A 

In the uplink, the subcarrier allocation and power allocation 
phases cannot be detached since subcarrier allocation depends on 
the powers constraints at every user. In Phase I, we define the term 

Xn,k,,l = -=te — Pn ' k i' 1 f or me a n oca tion of resources to 

users. This criterion guarantees the selection of the users who have 
not only better power-gain product but also result in less interference 
to the neighbour cells as the denominator J2j=i j^i fjsj.maxffn.fcj.U 
accounts for the maximum interference that the kf 1 user may cause to 
cell j. The details are provided in Algorithm 2. Once the subcarrier 
allocation is performed, the optimal powers can be calculated for 
the high SINR regime as explained above or for the general SINR 
regime through solving a series of GPs using successive convex 
approximation which is a provably convergent heuristic [121. This 
approach is known to compute globally optimal power allocations in 
many cases. Nevertheless, like any sub-optimal scheme, the conver- 
gence and optimality of the power allocation phase in the general 



SINR regime depends on the initialization of the power allocation 
variables. However, we observe that fast convergence can be achieved 
by using the optimal powers computed in the high SINR regime for 
initialization. The details of the implementation can be found in the 
extended journal version 1131 of this paper. 

Algorithm 2 Centralized Resource Allocation for L Cells 

1) Input: [H;], [A;], [Pi], [Gy] where a„, fei , ; = 0,p n , fei .i = 

[P kl ,mex/N] Vfc,,Vn 

Subcarrier Allocations (Phase I) 

2) Compute Xn,k L ,i for every k\ user at 71 th subcarrier in cell I. 
n = 1, do while n < N, n = n + 1 

3) Find the (n,ki,l) pair with maximal value of Xn,k t ,l- 

4) Perform allocations at n th subcarrier in the remaining cells 
jf, Vjr 7^ I by searching for the user with maximum value of 

Xn,k t ,j- 

5) Perform power update (equalization) as detailed in Algorithm 
1 in each cell I. 

6) Remove the subcarrier n from the set of unallocated subcarri- 
ers. end 

Power Allocations (Phase II) 

7) Compute the optimal powers P; in the high SINR regime using 
(10) given the allocations from Phase I. 

Complexity Analysis: Initially we perform a three dimensional 
search over Xn,k h j which has a complexity of O(LKN). Next, a 
linear search is conducted in the cells other than the starting cell 
with a complexity of 0((L — 1)K). The process iterates until all 
subcarriers are allocated, thus the total complexity of Phase I is 
0(LKN 2 + NK(L - 1)). The complexity of Phase II is difficult 
to determine, however, it can be measured in terms of degree of 
difficulty (DoD) that in turn relies on the number of constraints and 
variables associated with the GP 1141 . Since we are dealing with LK 
power constraints and LKN power variables, apparently it seems that 
implementing centralized GP/successive GP based schemes may not 
be a good choice for practical implementations. However, in order 
to reduce the complexity and DoD of GP we have developed the 
following less complex centralized scheme. 

B. Centralized Scheme B 

As the time complexity of Phase I is delectable compared to Phase 
II, in this approach, instead of splitting the whole procedure into 
two phases we combine them in the following manner. At first, the 
allocations are performed in each cell I followed by the power update 
(based on equalization) at subcarrier n as mentioned in Algorithm 
2 (Steps 2-5). Before performing the next subcarrier allocation, we 
compute GP based optimal power allocations on subcarrier n as the 
powers dictated by the equalization process Pn,k lt i eil may not be 
optimal. Setting the powers p n ,ki,l e as the upper bound on the power 
variable of fc* h user in each cell / at subcarrier n, we now define the 
less complex GP problem in order to compute the optimal powers 
that can maximize the throughput at n th subcarrier. 

L 



minimize 



io g2 r 



a 2 + I n 



^ - \ Pn,ki ,1 hn,ki ,1 

subject to p n ,k,,i < Pn,b,l„, VZ 



(11) 



Clearly, optimal power of each user at subcarrier n may not succeed 
in reaching the upper bound. Hence, this unused proportion of 
power that a user cannot use in reaching the upper bound, must 
be taken into account while doing the next round of allocations. 



Thus, this unused proportion of power is distributed among the 
remaining unallocated subcarriers of that user. The DoD reduces to 
L constraints and variables. Although this procedure restricts the 
degree of freedom offered by GP, numerical results show that the 
network throughput remains comparable with reduced complexity. 
The procedure is detailed in Algorithm 3. 

Algorithm 3 Low Complexity Centralized Resource Allocation for 
L Cells 

1) Compute Xn,k t ,i for every fc* h user at n th subcarrier in cell I. 
n — 1, do while n < N, n — n + I 

2) Steps 2-5 of Algorithm 2. 

3) Compute the optimal powers Pn,fc,,i in the high SINR regime 
using (11). 

4) Perform power equalization of kj h user among the remaining 
unallocated subcarriers. 

5) Remove the subcarrier n from the set of unallocated subcarri- 
ers. end 



V. Semi-distributed Resource Allocation Scheme 

We assume each BS knows the maximum interference caused by 
its users to the first tier of interfering cells. Assuming the knowledge 
of interfering gains here rendered this approach semi-distributed; 
however, we present in [13] a more distributed scheme that does 
not require the sharing of interfering gains among BSs. In the 
proposed semi-distributed scheme, each BS searches for the user 
with maximum value of Xn,ki,l m order to perform allocations. The 
allocation decisions are locally made at each BS (Phase I) and do 
not depend on the resource allocations and powers in the neighbour 
interfering cells. Once the allocations are decided each cell shares 
them with all other interfering cells. The optimal distributed power 
allocations can then be computed using dual decomposition methods 
by first performing the log transformation of the variables, i.e., 
In Pn,k t ,l = Pn"k,,i and In p n , kj ,j = Pn"k } ,j, then adding auxiliary 
variable z U} ij = Y^,j=i,j^i hn.kjjPn.kjd in order to transfer the 
coupling in the objective to coupling into the constraints 1121 . This 
reduces the complexity of implementing GP by decoupling the primal 
problem into L dual problems, i.e., one for each cell. The problem 
can be then be written as follows: 

L Ki N 

minimize ^ ^ ^ a n ,)t,,i log 2 



1 = 1 fc, = l n=l 

jV 



subject to } y a n ,hi,lPn,hi,l < Pfy.max, Vfc>,V7 (12) 

71=1 

L 

e <« > h n , kj ,ie Pn ' l i'i , Vn 

Writing the Lagrangian L(p„" kl ,i, z^ij, \ kl , Vn,ij) for (12) and split- 
ting the problem into L sub-problems yield (13) with local variables 
Pn,fcj,i, Zn,lj> Afc, and coupling variables T] n> ij- The details of the 
solution are omitted here due to space limitation but presented in 
the journal version of this paper 1131 . 

max min L(p n "k h t, Zi,lj> Vn.lj) (13) 

VI. Simulation Results and Performance Evaluation 

A cellular OFDMA network is considered where the radius of each 
cell is assumed to be R c = 0.5km. The users lie at equally spaced 
angles from 0 to 27r. For demonstration purpose, the total number 



of users are assumed to be same in all the cells. The maximum user 
transmit power is considered to be 2W. The channel gain is defined as 
h„ M ,i = (-122-107logi 0 rffc ! , i )-A / '(0,(7 2 ) + 101og 10 F„, fci , i . The 
first term denotes the path loss where 7 is the path loss exponent and 
is set equal to 3.8. The second term represents log-normal shadowing 
with a mean of 0 dB and a standard deviation of 8 dB. The last factor, 
F n ,ki,i corresponds to Rayleigh fading. The bandwidth of the system 
is assumed to be 5 MHz with a noise power spectral density of 2.07 x 
10" 14 W/Hz at each receiver. The channel conditions are assumed to 
be fixed during a frame. The interfering gains from j th interfering cell 
to the cell of interest I are computed as follows g n ,k-,i = ( — 122 — 
10 7 logi 0 d fcj , ; ) - Af(0, a 2 ) + 10 log 10 F„, fcjli 

In Table 1, we compare the maximum achievable average network 
throughput (measured in bps/Hz) and complexity of the proposed 
centralized and semi-distributed schemes with the derived bounds and 
the optimal solution in the high SINR regime. The optimal solution is 
computed via an exhaustive search based subcarrier allocation phase 
detailed in Section IV. The results are taken after averaging over 100 
channel realizations. It is observed that the performance gap of the 
centralized strategies and the optimal solution is delectable. However, 
this outcome may not remain valid for higher dimensions. In addition, 
the less complex centralized scheme B provides comparable results 
to centralized scheme A. It is also worth to mention here that in order 
to get a more tighter LB, we compute the LB throughput in Fig.l 
and Fig. 2 by computing the optimal power allocations based on GP 
given the LB allocations from Algorithm 1. Next, we compare the 
performance of the centralized scheme B and the semi-distributed 
scheme with the UB and LB for varying number of users and 
subcarriers in 2-cell (See 1131 ), 4-cell and 7-cell scenario. The results 
have been taken after averaging over 1000 channel realizations. We 
observe that for the 2-cell case, the centralized and semi-distributed 
strategies give nearly similar results, however, as the number of cell 
increases the performance gain of centralized scheme over semi- 
distributed scheme is evident (see Fig.l and Fig. 2). The UB results 
remain independent to the number of cells, however, it is clear from 
Fig.l and Fig. 2 that the LB throughput continues to decrease with 
the increase in number of interfering cells. Moreover, it can be easily 
observed that the performance gap between the centralized and semi- 
distributed strategies tends to increase with the increase in number of 
subcarriers, users and interfering cells. The established UB and LB 
provide a comprehensive idea of a finite performance gap. 

Another important observation in GP based optimal power al- 
location is that the sum Y^ n =iP n > k i> 1 ma y or ma y not ^ e eQua l 



to P, 



ki ,max 



as dictated by the equalization based power allocation. 



Moreover, it may not (almost never) be optimal that two users 
transmit with powers suggested by power equalization based power 
allocations. Consider a channel realization Hi,H2= [0.3028 0.2503; 
0.0388 0.1529]xl0~ 9 and Gi 2 ,G 2 i = [0.0597 0.0455; 0.1555 
0.0597]xl0~ n and [0.1340 0.6860; 0.7538 0.1935]xl0~ n , respec- 
tively. The equal power allocations dictates Pi = [0 0.5; 0 0.5] 
and P2 = [0.5 0; 0.5 0] for semi-distributed scheme that leads to 
the average throughput of 1 1.8392bps/Hz. However, computing the 
optimal powers results in Pi = [0.0000 0.5326; 0.0000 0.4674] 
and P 2 = [0.3840 0.0000; 0.6160 0.0000] which lead to a max- 
imum average network throughput of 17.2734bps/Hz. This example 
demonstrates the significance of GP as well as centralized scheme 
B in which after allocating each subcarrier, optimal powers are 
computed and thus, the remaining unused power can be accounted 
for while computing the next allocations. The significance of GP 
is also demonstrated in Fig. 3, where we compare the performance 
of the semi-distributed scheme with and without GP based power 



TABLE I 

Average network throughput of the centralized/semi-distributed schemes for L=2 cells and TV = 6 subcarriers/cell 





K=2 


K=4 


K=6 


Complexity 






di,fc,,!=0.35km 


di,fe,,!=0.45km 


d ii yt i ,i=0.35km 


*,fc„!=0-45km 


di,fc,,i=0.35km 


<2i,fc,,i=0.45km 






UB 


52.6264 


30.5791 


63.5723 


40.4159 


69.1927 


46.1368 


0(KN 2 ) 




Optimal 


39.3301 


20.0545 


49.6299 


28.4265 


55.6299 


34.8121 


0(K NL ) + 


DoD(LKN) 


Centralized A 


38.6454 


19.2522 


48.7848 


27.5679 


54.5671 


33.8055 


0(LKN 2 - 


- (L - \)KN) + DoD(LKN) 


Centralized B 


38.5254 


19.0222 


48.5248 


27.5479 


54.3671 


33.7925 


0(LKN 2 - 


- (L - \)KN) + DoD(L) 


Semi-distributed 


37.9408 


18.7190 


48.0243 


27.5107 


54.5178 


33.7270 


0(KN 2 ) + 
0(KN 2 ) 


DoD(KN) 


LB 


35.2340 


16.3582 


43.6006 


22.3586 


48.3627 


27.4565 





allocations. Without GP, the resulting throughput is observed to be 
lowered significantly than the centralized scheme. 

VII. Conclusion 

In this paper, we developed an UB and LB to the problem 
of sum rate maximization in multi-cell uplink OFDMA networks. 
Moreover, we proposed two centralized and one semi-distributed 
resource allocation schemes. All schemes are evaluated and compared 
in terms of network throughput and computational complexity. 
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Fig. 1. Comparison of the centralized B/semi-distributed schemes for L—l 
cell, d=0.4km 
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Fig. 2. Comparison of the centralized B/semi-distributed schemes for L—A 
cell, d=0.4km 
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Fig. 3. Significance of GP based optimal power allocations in high SINR 
regime for L=2 cells, 7V=10 subcarriers/cell and d=0.4km: An example 



